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The survival probability, Sw(t), of a diffusing prey ("lamb") in the prox- 
imity of N diffusing predators (a "pride of lions") in one dimension is 
investigated. When the lions are all to one side of the lamb, the survival 
probability decays as a non-universal power law, Sjv(t) oc t~@ N , with 
the decay exponent 0n proportional to IniV. The crossover behavior 
as a function of the relative diffusivities of the lions and the lamb is 
also discussed. When N —>■ oo, the lamb survival probability exhibits a 
log- normal decay, Soo(t) cx exp(— In 2 t). 



I. INTRODUCTION 



Consider a particle system which consists of a diffusing prey and iV independent, 
diffusing predators, with iV either finite or infinite. The prey is absorbed, or dies, when- 
ever it is touched by any of the predators. We are interested in the probability Sjv(£) 
for this "lamb" to survive until time t when it is besieged by these predatory "lions" [1]. 
While this appears to be a simple problem, there are important aspects of the long-time 
behavior which are incompletely understood. Their resolution has fundamental ramifi- 
cations for diffusive processes in the presence of complex absorbing boundaries and also 
practical implications, as this type of capture process appears in a variety of applica- 
tions, such as diffusion-controlled chemical kinetics, wetting, melting, and commensurate- 
incommensurate transitions (see, e. g., Refs. [2-4]). It is known rigorously that for spatial 
dimension d > 3, the capture process is "unsuccessful" (in the terminology of Ref. 1), in 
that there is a finite probability for the lamb the survive indefinitely for any N and for any 
initial spatial distribution of the lions. This result is a consequence of the transience of ran- 
dom walks for d > 2 [5]. For d = 2, the capture process is "successful" - the lamb dies with 
probability one. However, diffusing lions in two dimensions are still sufficiently poor preda- 
tors that the average lifetime of the lamb is infinite. Moreover, Sjv(t) oc [Si(t)] oc [lnt] _7V , 
i. e., the many-body nature of the capture process is basically irrelevant. 

In one dimension, however, diffusing lions are more efficient in their predation because 
of the recurrence of random walks [5]. This typically leads to a lamb survival probability 
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which decays as a power law in time. More generally, we may consider the survival prob- 
ability, SN ; M{t) = SM,N(t), when initially N lions are placed to the left and M lions are 
placed to the right of the lamb. The most interesting situation is when initially N lions are 
to one side of the lamb. For this initial condition, realizations in which the lamb runs away 
from the lions leads to an anomalously slow decay of the lamb survival probability. Our 
primary result is to determine the asymptotic behavior of the lamb survival probability in 
this predatory environment, Sjv,o(£) = <SW(£)- For finite N, we predict that SN(t) ~ t~^ N 
with /3n oc In N, while 5' 00 (t) oc exp(— In 2 t) for N — > oo. We shall argue that these depen- 
dences arise from the fact that the motion of "closest" lion (whose individual identity can 
change with time) is enhanced compared to normal diffusion. When the number of lions is 
finite, this enhancement manifests itself in the diffusivity of the last lion being proportional 
to IniV. This factor is ultimately responsible for the logarithmic increase of [3n on N. In 
the limit N — > oo, the co-ordinate of the last lion actually varies as y/thit, and leads to 
the lamb survival probability varying as exp(— ln 2 t). 

To provide context for our results, consider first a system consisting of one lamb and 
one lion, i. e., (M,N) = (1,0) or (0, 1). The survival probability is trivially calculable in 
this case, since the distance between the lion and the lamb undergoes pure diffusion with 
an associated diffusivity Dl + Dg. Here Dl and Dg are the lion and lamb diffusivities, 
respectively. Because of the equivalence to diffusion, the survival probability is [5] 

S 10 (t)~^= X ° (1) 
^{D L +Dg)t 

where xq is the initial separation between the lion and the lamb. A more interesting situa- 
tion is that of two lions with either: (i) one lion on either side of the lamb ( "trapped" lamb) , 
or (ii) both lions to one side ("chased" lamb). These two systems can be straightforwardly 
solved when the diffusivities of all three particles are different. For the trapped lamb, 
denote the particle positions as x±, X2, and x 3 , with 1 and 3 referring to the surrounding 
lions, and 2 to the trapped lamb. Let the corresponding diffusivities be D±, D 2 , and D3. 
To solve for the survival probability, it is convenient to introduce the rescaled co-ordinates 
Hi = Xi/VD~i, each of which diffuses at the same rate. The survival of the lamb corre- 
sponds to the constraints y\\fD\ < y 2 \[D 2 ~ and y 2 \/D 2 ~ < yz\/Dz- Since the co-ordinates 
yi diffuse isotropically, lamb survival is equivalent to the survival of a random walk in 
three-dimensional space within the wedge-shaped region bounded by the absorbing planes 
ywfiJ\ = y 2 \[D 2 and = 2/3 v^3- By straightforward geometric considerations [2], 

this three-dimensional diffusion process is equivalent to diffusion in two dimensions within 
a wedge of opening angle 

9 = cos" 1 ( ° 2 =) (2) 

\ V /(D 1 +D 2 )(D 2 +D 3 ) J 

For this two-dimensional problem, it is well-known that the survival probability asymp- 
totically varies as t -71 "/ 26 * [6]. Identifying Di = D 3 = D L and D 2 = Dg, leads to 
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where TZ = D^/D L . Similarly, one finds for the "chased" lamb 

(36) 

Physically, when TZ — > oo, the motion of the lions becomes irrelevant and the trapped 
lamb problem maps to a diffusing particle in a fixed size absorbing domain. For this 
geometry, the survival probability decays exponentially in time, corresponding to — > oo 
in Eq. (3a). Conversely, for the lamb at rest, TZ — > 0, the survival probability is the square 
of the corresponding survival probability in the two-particle system so that /^i — > 1. 
For the chased lamb (Eq. (3b)), /?2,o has the limiting values /?2,o = 1/2 for TZ — > 00 and 
/?2,o = 1 for 7?. — > 0. These two values are in accord with a direct consideration of the 
extreme cases of immobile lions or immobile lamb, respectively. Notice also that for TZ = 1 
we get = 3/2 and /?2,o = 3/4, results which were known previously [1]. 

Extending the above approach, the survival of a lamb in the presence of N > 2 lions 
can be mapped onto the survival of an N + 1-dimensional random walk which is confined 
within an absorbing hyper-wedge. This problem does not appear to be generally soluble, 
however. Numerically, there has been investigation [1] of the one-sided equal-diffusivities 
problem for the cases iV = 3, 4, and 10 as part of an effort to understand the general 
behavior on N. This simulation reveals that the exponent /3n grows slowly with N, with 
Ps = 0.91, (3 4 ps 1.032 ± 0.01, and (3 10 ps 1.4. (Because the case N = 4 is close to 
the transition between a finite and infinite lamb lifetime, there is more numerical data 
and hence a greater precision in the estimate for (3^.) The understanding of this slow 
dependence of (3n on iV is the focus of our work. 

In the next section, we provide a heuristic argument for the dependence of (3n on N, 
as well the behavior for N — > 00. A more complete derivation of these results is given Sec. 
III. The general dependence of the survival probability on the ratio TZ = D^/Dl is also 
considered. In Sec. IV, we treat the case of N — > 00. A general discussion and conclusions 
are given in Sec. V. 

II. HEURISTIC ARGUMENTS FOR ONE-SIDED SIEGE 



2 cos" 1 

7T 1 + 1Z 



First consider the trivial case of a stationary lamb, Di = 0. For non-interacting lions, 
the lamb survival probability is just the product of the survival probabilities associated with 
each lion-lamb pair. This immediately gives Sjv(£) oc (t~ l / 2 ) N , from which /?jv(0) = N/2. 
For this case, the relative positions of the lions do not matter in the asymptotic behavior 
of S'jv(t), i. e., it is immaterial whether the lamb is "trapped" or "chased". When the lamb 
also diffuses, it is convenient to consider the survival probability in the rest frame of the 
lamb. Although the lions still diffuse independently, their relative motions with respect to 
the lamb are not independent. Therefore to determine the survival probability of a diffusing 
lamb, it more useful to track the position of the closest lion only. For concreteness and 
simplicity, suppose initially that all the lions are at the origin and the lamb is at xq > 0. 
A rough estimate for the location of the closest or "last" lion, is provided by 

1 -1 

- e -* 2 /4D L t dx = 1/Ar _ (4) 

x+{t) V4nD L t 
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This specifies that there should be one lion in the range (x + (t), oo) out of an initial group 
of N lions. By asymptotic expansion of this integral, the location of the last lion is given 
by 

x+(t) ~ y/4D L t InN. (5a) 

In the limit N — > oo, a physically tractable initial condition is to have a uniform density 
of lions Co extending from — oo to 0. In this situation, only a number N oc yj c^D^t of 
the lions are "dangerous", i. e., are potential candidates for being closest to the lamb. 
Consequently, for N — > oo, the leading behavior of x + (t) is 



x+(t) ~ \j2D L t \n(clD L t) N -> oo. (56) 

The next step in our heuristic approach is to posit that for large N, the true stochastic 
motion of the last lion can be replaced by a continuous motion x(t) with x(t) = x+(t), 
as given in Eqs. (5). Then the system reduces to a two-body problem of a diffusing lamb 
and an approaching absorbing boundary, whose location is x + (t) = y/At. As discussed 
in the next section, the survival probability of a lamb adjacent to such an approaching 
"cliff" can be calculated analytically. This gives the exponent of the survival probability 
as (3 ~ A/16Dg. Substituting the appropriate value of A as specified by Eqs. (5), we obtain 
S N (t) ~ t-/M^) w ith N (U) ~ ]n(NK)/4fc for finite N, and S^t) ~ exp(-ln 2 t) for 
N -> oo. 

III. ASYMPTOTIC ANALYSIS FOR ONE-SIDED SIEGE 

A more rigorous approach is to consider the survival probability in an auxiliary "dead- 
line" problem whose asymptotic behavior turns out to give a tight lower bound for the 
true survival probability of the lamb. The deadline problem is defined as follows: Consider 
an imaginary point Xdead(^) between the lamb and the lions which moves deterministically 
according to Xd ea d(t) = y/At. If the lamb crosses this line, it is considered to have died; 
analogously, if any of the lions overtakes the deadline, the lamb is again considered to 
have died. Our strategy is to determine the survival probability in this auxiliary problem, 
and then maximize this probability with respect to the parameter A. First notice that a 
deadline position which is proportional to y/i already optimizes the lamb survival proba- 
bility with respect to other power-law motions for the deadline. That is, is xa ea d(t) were 
proportional to t a with a < 1/2, we would asymptotically recover the behavior for the 
stationary deadline, which grossly overestimates the decay exponent as [3n = (N + l)/2. 
Conversely, for faster than diffusive deadline motion, i. e., a > 1/2, the probability that 
the lamb does not hit the deadline decays as a stretched exponential [7]; therefore, this 
case can be also ignored. The marginal situation of a = 1/2 thus plays a fundamental role. 

To compute the survival probability for the deadline problem, we have to solve two 
first-passage problems: (i) The survival of a diffusing particle in the proximity of a receding 
absorbing boundary, or cliff. This corresponds to a single lion, and we define the probability 
that a lion does not reach the cliff to be Sn on (t) oc t~^ lion ^ Dl " A \ (ii) The survival of a 
diffusing particle in the proximity of an approaching cliff. This corresponds to the lamb, 
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and the associated survival probability is defined as S\ am b(t) oc t~^ aiah ^ De ' A \ rp^g f u j[ 
survival probability S^it) oc t~^ N is clearly the product S^ OQ (t) 5i am b(t), so that (3n = 
Aamb(-D^ 5 A) + N /3iion(-DL, ^4) • Once we know the exponents (3u on and Aamb, we optimize 
the decay exponent (3n with respect to the amplitude A. By appealing to the method of 
the "optimal fluctuation" [8] , we hypothesize that this extremal survival probability in the 
deadline problem gives the true asymptotic behavior. 

Fortunately, the exponents Aamb and /3ii 0n have been computed in various physical 
[9] and mathematical studies [fO] so that the full deadline problem is soluble. For com- 
pleteness, however, we outline our approach, given in Ref. [7], which has the advantage of 
conceptual and technical simplicity. While this earlier work considered only the case of 
a receding cliff (relevant for the lions), the extension to the case of an approaching cliff 
can be derived with minimal additional effort. Let us therefore recall the steps in the 
computation of the survival probability for the case of the receding cliff. Consider a lion 
which is initially placed on the negative x-axis and that the cliff position is xo(t) = \f~At. 
In the long-time limit, the lion density approaches the scaling form [7] 

c(x,t) ^r^-^cco, (6) 

where £ = 1 — — is the appropriate dimensionless distance variable and C(£) is a scaling 
function. The initial co-ordinate of the lion, — oo < x < xq, corresponds to < £ < oo. 
The power law prefactor is chosen to ensure that the survival probability decays as t~™ on , 
as defined previously. 

Substituting Eq. (6) into the diffusion equation, one finds that the scaling function 
satisfies 

D L d 2 C 1,. ^dC ( 1\. n 

+ o(C-l)^+ Aion + 77 C = 0. (7) 



A di 2 2^ ' di 
Introducing the transformation 



Z- 1 = \l^r^ C(0=exp(-^j2%), (8) 
one finds that T>(rj) satisfies the parabolic cylinder equation of order 2(3u on [12], 



+ 



2Aion + 2 ~ 4 



^2/3 lion = 0. (9) 



dr] 2 

The absorbing boundary condition at the edge of the cliff implies 

V (-y/A/2D L } = 0. (10a) 

On the other hand, to avoid a singular solution at r\ = oo, the second boundary condition 
is 

V(rj = oo) = 0. (106) 

Mathematically, the determination of (3u on and T>(rj) is equivalent to finding the ground 
state energy and wave function of a quantum particle in a potential composed of an infinite 
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barrier at r\ = — ^JA/2D L and the harmonic oscillator potential for r\ > — ^JA/2D L 
[11]. Higher excited states do not contribute in the long time limit. This relation with 
quantum mechanics allows one to apply well-known techniques to determine the asymptotic 
behavior [7,10]. Among the two elemental solutions of the parabolic cylinder equation, 
^ > 2/3 lion (??) and V 2/3 lion (— 77), only the former satisfies the boundary condition V{oo) = 0. 
Therefore, the absorbing boundary condition of Eq. (10a) determines the decay exponent 

Aion = Aion(-D,L, A). 

As discussed previously, the interesting behavior emerges in the large N limit. For 
this case, the deadline position grows as \J~A~i but with an anomalously large amplitude A. 
Consequently, the probability distribution of each lion is only weakly affected by the reced- 
ing deadline. This allows us to employ the "free particle" particle Gaussian approximation 
for the probability distribution of each lion. Although this form does not satisfy the ab- 
sorbing boundary condition, the error is negligible because A>1. Consequently, we can 
determine the decay exponent simply by computing the flux to the absorbing boundary 
for the assumed Gaussian probability distribution [7]. This yields 



In the limit A — > 00, this simple-minded approach coincides with the results from a com- 
plete analysis in terms of the parabolic cylinder function solution. 

An analogous, but simpler, treatment applies for the approaching cliff, which we use 
to describe the interaction of the deadline with the lamb. That is, suppose that a lamb is 
initially placed on the positive x-axis and that there is an approaching cliff whose location 
is at \^At. To solve this problem by the same approach as the receding cliff, we introduce 
the appropriate dimensionless length variable £ = — 1 and make the analogous scaling 
ansatz as in Eq. (6), so that Eq. (7) is replaced by 

T? + i« + i)f + (/wi)c = o. (12) 

For this case, it is again helpful to introduce rj via £ + 1 = \J ^j 1 rj and C(£) = e~ v / 4 T>(rj). 

The scaling function T>(rj) is again the parabolic cylinder function of order 2f3\ am h and the 
absorbing boundary condition, 

V2^ h (y/A/2D e ) =0, (13) 

now determines the decay exponent (3\ am b = f3\ am b(Dg, A). Since the relevant zero of the 

parabolic cylinder function 77 = a/ A/2Dg is large, Aamb is also large. Then an inspection 

2 

of Eq. (9) provides the estimate 2/9i am b + \ ^ \, or 

p lamh (D e ,A)~-£-. (14) 
Therefore the total decay exponent for the deadline problem is 



p N (K, A) = Aamb(^, A) + N[3 lion (D L , A) ~ ^A- + N^^- e~ A l^. (15) 
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Minimizing this expression with respect to A yields the optimal value A* ~ AD L \n(AN1Z). 
Thus the deadline motion is enhanced by a factor of In N compared simple diffusion; notice 
that this coincides with the motion of the last lion in a pride of N lions. Correspondingly, 
the optimal value of the decay exponent (3n(JZ) = Pn(TZ, A*) is 

MK) „ (16) 

Our construction of the deadline problem relies on the assumption that N >> 1. This 
assumption is crucial, otherwise the deadline problem would not provide a meaningful 
approximation for the behavior of the original system. However, the physical nature of 
the problem suggests that different asymptotic behaviors for the lamb survival probability 
should arise for 1Z ^> 1 and 7£ <C 1. In fact, consideration of the limiting cases of a 
stationary lamb and of stationary lions, suggests that Eq. (16) is actually valid only for 
N~ x <C 1Z <C InN. In the slow-lamb limit, TZ <C iV -1 , the logarithmic behavior of 
Eq. (16) should cross over to that of the stationary-lamb limit, namely, /3/v(0) = N/2. In 
the complementary fast-lamb limit, TZ ^> InN, the behavior of the stationary lion case 
should be recovered, in which /?jv(oo) = 1/2. Thus the full dependence of (3n on the 
diffusivity ratio 1Z is expected to be 

IN/2 ft < 1/iV; 

\n{ANK)/AlZ l/N <^TZ < IniV; (17) 

1/2 ft>lniV. 

The non- universal dependence of (3n on the diffusivity ratio for the intermediate regime of 
l/N <C TZ <C InN is the generalization of the exponents in Eqs. (3), for the three-particle 
system, to arbitrary N. 

IV. INFINITE NUMBER OF LIONS 



Consider now a lamb which is under one-sided siege by an infinite pride of lions. (These 
lions need to be distributed over an infinite domain so that their density is everywhere 
finite. If the lion density were infinite at some point, then the closest lion would move 
inexorably toward to the lamb at each step, leading to the survival probability decaying 
exponentially in time.) The interesting situation is when the lions are all to one side of 
the lamb. However, to introduce our approach, it is instructive to consider first the simpler 
two-sided problem, in which lions are uniformly and symmetrically distributed with unit 
density on either side of a stationary lamb, a problem has been previously investigated by 
asymptotic and exact methods [13,14]. For completeness, we describe an approach which 
is in the spirit of the previous section. 

For a lamb at the origin, the density of the lions c(x, t) may be found by solving the 
diffusion equation with an absorbing boundary condition at x = and with the initial 
condition of a unit density everywhere. This yields [5] 
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2 px\/vmzi _ 2 
ciamb(M) = —= I d(e~ c 
V n Jo 



dc 



Thus the diffusive flux of lions toward the lamb is -Dl(^| x=0+ 
survival probability 5"^ (t) therefore obeys 



dc 

dx i a 



.0-) = y 



4Dl 

TTt 



dSUt} - -*»(«) IADl 



with solution 



eft 



S'ooW = exp 



nt ' 



-4 



7T 



(18) 



The 



(19) 



(20) 



When both the lion and lamb are diffusing, a faster decay occurs. However, since the 
dominant annihilation mechanism arises from the diffusive flux of lions toward the lamb, 
we expect that the asymptotic decay is still given by Eq. (20) [13]. The crucial feature 
of this two-sided problem is that there is no good "survival" strategy, so that the lamb 
survival probability must decay rapidly in time. 

For a lamb under one-sided siege, we again attempt a solution via the auxiliary 
deadline model. Assume that the deadline undergoes enhanced square- root motion, i. e., 
Xdead(t) = y/~At, with A ^> 1. Repeating the steps employed previously for a finite pride of 
lions, we have S' 00 (t) = >S'iion(£)<S , iamb(£), with S\ am h(t) oc t~ A / 16De , as in the case of a finite 
pride. To determine Sy lon (t) we again use a free particle approximation, since we expect 
that the amplitude A will be large. Thus for the probability density of the lions, we ignore 
the adsorbing boundary condition on the moving deadline. For the initial condition of unit 
density of lions for x < and zero density otherwise, the time dependent lion density is [5] 



Clion(^-? t) 



1 



d(e 



-C 2 



(21) 



Although this solution disagrees with the adsorbing boundary condition on the deadline, 
the disagreement is of order e~ A / 4DL and is negligible when A ^> 1. 

Computing the diffusive flux of lions through the deadline, we make use of Eq. (21) 
and A/D L > 1 to find 



D ^° 
L dx 



X — 2^dead (^) 



Dl ( 



-A/AD! 



(22) 



The lion survival probability, Su on (t), therefore obeys 

dSy lon (t) 



dt 



with solution 



5iion(*) - exp 



-A/ADl 



Dr.t 



7V 



(23) 



(24) 
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Thus the full survival probability is 



(t) oc exp 



-A/AD L 



Drt 



7V 



A 



Int 



(25) 



Maximizing this survival probability with respect to A, we find that optimal value, 
A* , grows in time as, 

A* ~ 2D L In 



ln 2 t 



(26) 



The leading logarithmic behavior is in accord with our naive estimate given in Sec. II. 
Combining Eqs. (25) and (26) gives 



Soo (t) oc exp [— In 2 t\ 



(27) 



Thus we obtain a survival probability for the one-sided system which decays faster than 
any power law and slower than any stretched exponential. The decay is, however, universal 
in that the power of logarithm does not depend on the diffusivity ratio TZ = D^/D^. 



V. SUMMARY AND DISCUSSION 



For a diffusing lamb in one dimension which is adjacent to a pride of iV diffusing, 
predatory lions, the survival probability of the lamb decays as S^it) ~ t~* N with (3n 
proportional to IniV. This slow increase of /3/v on N reflects the fact that the dominant 
contribution to the survival probability arises from realizations in which the lamb "runs 
away" from the lions. Consequently, each additional lion in the system has a progres- 
sively weaker effect on the lamb survival. This is in contrast to the case of a stationary 
lamb, where each additional lion is equally effective in hunting the lamb, so that (3n is 
proportional to N. The exponent @n is also a decreasing function of the diffusivity ratio, 
K = D £ /D L , with p N = N/2 for K = and (3 N = 1/2 for K = oo. Thus, in accord 
with intuition, the best survival strategy for the lamb is to diffuse faster than the lions. 
In contrast, for a two-sided system, where the lions initially surround the lamb, the best 
survival strategy for the lamb is to remain still. 

The above non- universal power-law decay of Sj^it) motivated the basic question, con- 
sidered in Ref. 1, of whether the mean lamb lifetime 

r N = - [°° dtt = r dt S N (t) (28) 

Jo dt J 

is finite or infinite. From Eq. (28), it is clear that r/v is finite for (3n > 1, and tn diverges 
otherwise. The numerical evidence from Ref. 1 indicates that when TZ = 1 the lamb lifetime 
is finite for N > 4. Since our prediction for [3n is anticipated to be accurate only for large 
N, we may conclude that the lamb lifetime is finite when iV > iV* (1Z) , but cannot provide 
an accurate estimate of this threshold value. Additionally, we predict that iV* (1Z) should 
increase rapidly with the diffusivity ratio Dg/D L: namely In N* (71) oc 1Z (Eq. (17)). 



9 



For an infinite number of predators, the lamb survival probability 5' 00 (t) exhibits a 
log-normal decay exp(— ln 2 t). This contrasts sharply with the corresponding behavior 
in the two-sided geometry, where 5' 00 (t) oc exp(— t 1 / 2 ). For the one-sided geometry, it is 
striking that this same survival probability occurs for a reactive system consisting of a 
single "fast impurity" which moves with velocity v > 1 within a semi-infinite sea of mutu- 
ally annihilating ballistic particles moving at velocity v = ±1 [15]. Given the superficial 
similarity of the fast impurity with the lions and lamb systems, it may be interesting to 
seek a deeper connection between these two problems. 

We close with mention of a generalization where the lions are "vicious" among them- 
selves, in addition to stalking the lamb. (For such self-predatory lions, their number 
must be infinite; otherwise, the lamb survival probability has a non-zero asymptotic 
value.) There are two natural possibilities for the outcome when two lions meet: either (i) 
one lion dies (aggregation), or (ii) both die (annihilation). The first possibility is particu- 
larly simple, since the closest lion undergoes pure diffusion, independent of its individual 
identity. Thus the two-sided geometry reduces to the finite particle system (M, N) = (1, 1), 
with decay exponent given by Eq. (3a). The one-sided aggregation problem is even simpler 
since it reduces to the (1,0) problem whose solution is given in Eq. (1). 

Annihilating lions framework leads to more interesting behavior, as the position of the 
closest lion suddenly jumps away from the lamb whenever the closest two lions annihilate. 
The two-sided version of this problem was introduced in Ref. 16. It was found that the 
lamb survival probability decays as a non-universal power law, S(t) oc t~ l( ^\ with a 
Smoluchowski theory predicting 7(72.) = a/(1 + 72)/8. This agrees with the obvious exact 
result 7(1) = 1/2 and is close to the exact value 7(0) = 3/8 [17]. This Smoluchowski 
prediction and also provides a good approximation for the simulation results for arbitrary 
72 [16]. Generalizations of the two-sided annihilation problems (e. g., to many dimensions) 
have also been discussed in [16-18]. To the best of our knowledge, however, lamb survival 
in the presence of a one-sided distribution of annihilating lions has not yet been treated. 
If one naively assumes that the two-sided death probability can be expressed in terms of 
independent one-sided death probabilities, then the exponents of the one-sided system, 
(3(H), and the two-sided system, 7(72), are simply related by (3(71) = 7(72)/2. This is 
clearly correct for 72 = 0, where the independence of the one-sided death probabilities 
is exact. Consequently, the known value of 7(0) gives (3(0) = 3/16. However for 72 > 
0, the motions of the lions are not independent in the rest frame of the lamb, and the 
independence of one-sided killing probabilities is only an approximation. Interestingly, 
however, simulations suggest that for equal lion and lamb diffusivities, (3(71 = 1) = 1/4, 
which conforms to the relation (3(71) = 7(72)/2. We do not have an understanding of this 
simple yet paradoxical result. Finally, as 72 — > 00, it is clear that (3 (TV) — > 1/2. Thus we 
conclude that (3(TV) is a slowly increasing function of 72, with (3(0) = 3/16 and (3 (00) = 1/2. 
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